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AN INDEX THEOREM FOR LORENTZIAN MANIFOLDS WITH 
COMPACT SPACELIKE CAUCHY BOUNDARY 

CHRISTIAN BAR AND ALEXANDER STROHMAIER 


Abstract. We show that the Dirac operator on a compact globally hyperbolic Lorentzian 
spacetime with spacelike Cauchy boundary is a Fredholm operator if appropriate boundary 
conditions are imposed. We prove that the index of this operator is given by the same 
expression as in the index formula of Atiyah-Patodi-Singer for Riemannian manifolds with 
boundary. The index is also shown to equal that of a certain operator constructed from 
the evolution operator and a spectral projection on the boundary. In case the metric is of 
product type near the boundary a Feynman parametrix is constructed. 


Introduction 

The Atiyah-Singer index theorem is one of main mathematical achievements of the 20**^ 
century due to its many applications and because of the conceptual insights it provides. 
It computes the index of an elliptic operator on a compact manifold without boundary. An 
analog for hyperbolic operators on Lorentzian manifolds is unknown and is not to be expected. 
Compact Lorentzian manifolds without boundary are for many reasons unsuitable in this 
context. For example, they violate all causality conditions, which not only makes them 
unsuitable as models in General Relativity but also causes analytic problems for the operators; 
e.g. there is no well-posed initial value problem. 

The Atiyah-Patodi-Singer index theorem deals with elliptic operators on compact mani¬ 
folds with boundary. The boundary conditions which one has to impose are based on the 
spectral decomposition for the induced operator on the boundary. The aim of the present 
article is to derive a Lorentzian analog for this theorem. 

We consider Lorentzian manifolds M with boundary where the boundary consists of two 
smooth spacelike Cauchy hyper surfaces. Since the boundary is Riemannian, Atiyah-Patodi- 
Singer boundary conditions still make sense. Moreover, we assume a spin structure on the 
manifold so that spinors and the Dirac operator are defined. In even dimensions the spinor 
bundle splits into the subbundles of spinors of positive and negative chirality, respectively, 
SM = S~^M(BS~M. The Dirac operator maps sections of positive-chirality spinors to those of 
negative chirality. The domain FE^pg{M; S~^M) of the Dirac operator can be characterized 
as the completion of smooth sections of S~^M with respect to the L^-graph norm for the 
Dirac operator, subject to Atiyah-Patodi-Singer boundary conditions. Theorems 3.3 and 4.1 
combine to give the 

Main Theorem. Let {M,g) be a compact time-oriented globally hyperbolic Lorentzian man¬ 
ifold with boundary dM = S_ U S_|_. Here are smooth spacelike Cauchy hypersurfaces, 
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with S_|_ lying in the future o/S_. Assume that M is even dimensional and comes equipped 
with a spin structure. 

Then the Dirae operator -Daps • S~M) under Atiyah-Patodi- 

Singer boundary conditions is Fredholm and its index is given by 

ind[DAPsl = f A(V) + / TA(9) - '■(A-) + ft(A,) + .,(A_)-,(A,) 

JM JdM ^ 

The right hand side in the index formula is precisely the same as in the original Riemann- 
ian Atiyah-Patodi-Singer index theorem. Here A(V) is the A-form manufactured from the 
curvature of the Levi-Civita connection of the Lorentzian manifold, TA( 5 f) is the correspond¬ 
ing trangression form which also depends on the second fundamental form of the boundary 
(and vanishes if the boundary is totally geodesic). Moreover, A± denotes the Dirac operator 
on h the dimension of the kernel, and r] the r/-invariant. 

The Dirac operator on a Lorentzian manifold is far from hypoelliptic; solutions of the Dirac 
equation can have very low regularity. Theorem 3.5 tells us that under Atiyah-Patodi-Singer 
boundary conditions this is no longer so. Solutions are always smooth as if we had elliptic 
regularity at our disposal. Moreover, 

ind[DAPs] = dimker[DAPs|c°°(Af;S+Af)] “ dimker[DaAPslc'°“(M;S+Af)] ■ 

Here DaAPS stands for the Dirac operator subject to “anti-Atiyah-Patodi-Singer” boundary 
conditions, the conditions complementary to the APS conditions. The occurrence of the 
aAPS boundary conditions and the fact that DaAPS again maps sections of S~^M to those 
of S~M, and not in the reverse direction, are different from the corresponding formula in 
the Riemannian setting. In the elliptic case, the Dirac operator with anti-APS boundary 
conditions will in general have infinite-dimensional kernel. In the Lorentzian situation, anti- 
APS conditions work equally well as the APS conditions. 

The Atiyah-Singer index theorem and the Atiyah-Patodi-Singer index theorem have often 
been used in Quantum Field Theory, particularly in the context of the chiral anomaly. Here 
the underlying space is assumed to be Riemannian and the formulae are then applied to 
Lorentzian spacetimes by analogy. This is partially motivated by perturbative computations 
of gravitational corrections to the chiral anomaly that yield the same terms as in the Atiyah- 
Singer index formula. Our theorem gives a direct and non-perturbative interpretation of the 
chiral anomaly together with an r/-correction term. 

The index of Dirac operators on curved Lorentzian spacetimes seems not to have been 
discussed much in the physics or mathematics literature. We believe that this is due to the 
lack of ellipticity of the operator and the fact that in the hyperbolic case regularity is implied 
by propagation of singularities of the boundary data rather than the existence of a local 
parametrix. We are aware only of [13], where an index is associated to a certain bounded 
operator acting on the solution space of the Dirac operator on some globally hyperbolic 
spacetimes. There seems to be no relation to APS-boundary conditions and to the results of 
the present paper. 

The paper is organized as follows: In Section 1 we describe the setup and collect some 
standard material on Lorentzian geometry and on Dirac operators, mostly to fix notation. 
In Section 2 we introduce spinors of finite energy and discuss well-posedness of the Cauchy 
problem in those function spaces. This gives rise to the wave evolution operator Q for the 
Dirac equation. The Atiyah-Patodi-Singer boundary conditions imply a decomposition of this 
evolution operator into a 2 x 2-matrix. We show that the off-diagonal terms Qh _and Q _ 
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are compact operators while the diagonal terms (5++ and Q _are Fredholm. This is based 

on a Fourier integral operator technique. Some technical material needed here is collected 

in Appendix A. In Section 3 we relate the Fredholm property of Q _to that of the Dirac 

operator itself and show that the indices coincide. Here we also show smoothness of solutions 
of the Dirac equation under APS or anti-APS boundary conditions. The geometric formula 
for the index is computed in Section 4. The spectral flow of a family of Dirac operators on 
Cauchy hypersurfaces is used to relate the index of the Lorentzian Dirac operator to that of 
a Dirac operator for an auxiliary Riemannian metric. In Section 5 we construct Lorentzian 
metrics on M = / x with nontrivial index. 

In Section 6 we describe an alternative approach to index theory for Dirac operators on 
Lorentzian manifolds based on the construction of a Feynman parametrix. We show that the 
Dirac operator subject to APS boundary conditions is invertible up to smoothing operators, 
at least in the case when the metric of M is of product type near the boundary. This shows 
that the theory in some ways still resembles the elliptic theory, but Fourier integral operator 
parametrices replace the pseudodifferential parametrices. This approach also shows that the 

operators _and Q _are smoothing (and not just compact) when the metric of M is of 

product type near the boundary. This is important for physical applications as it implies the 
implement ability of time evolution in the Fock space constructed from the space of solutions 
of the Dirac equation, by the Shale-Stinespring criterion. 

We conclude in Section 7 with some remarks on possible extensions of the results of this 
paper. 

Acknowledgments. The authors would like to thank the Sonderforschungsbereich 647 
funded by Deutsche Forschungsgemeinschaft, and the London Mathematical Society (Scheme 
4 grant 41325), for financial support and the universities of Loughborough and Potsdam for 
their hospitality. Part of the work for this paper was carried out during the trimester pro¬ 
gram “Non-commutative Geometry and its Applications”, funded by the Hausdorff Center of 
Mathematics in Bonn, and we are grateful for the support and hospitality during this time. 

1. Setup and notation 

We start by describing the setup of this article and collect a few standard facts on Dirac 
operators on Lorentzian manifolds. For a more detailed introduction to Lorentzian geometry 
see e.g. [8, 18], for Dirac operators on semi-Riemannian manifolds see [5, 6]. Suppose that 
X is an (n -|- l)-dimensional oriented time-oriented Lorentzian spin manifold. We use the 
convention that the metric of X has signature (—h ■ ■ ■ -|-). 

A subset S C A is called a Cauchy hypersurface if every inextensible timelike curve in 
X meets S exactly once. If X possesses a Cauchy hypersurface then X is called globally 
hyperbolic. All Cauchy hypersurfaces of X are homeomorphic. We assume that X is spatially 
compact, i.e. the Cauchy hypersurfaces of X are compact. Let C M be two disjoint 

smooth and spacelike Cauchy hypersurfaces. W.l.o.g. let S_ lie in the past of By [9, 
Thm. 1.2]^ X can be written as 

(1) A = ]RxS 


^The cited theorem ensures that one given Cauchy hypersurface can be prescribed as a level set of the 
Cauchy temporal function. But its proof is local and can easily be adapted to the case of two disjoint Cauchy 
hypersurfaces, here E±. 
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such that each = {t} x S is a smooth spacelike Cauchy hypersurface, S_ = and 
= Etj. The metric of X then takes the form (•, •) = —N'^dt^ + gt where : X —)• M 
is a smooth positive function (the lapse function) and gt is a smooth 1-parameter family 
of Riemannian metrics on S. We understand the subset M = [ti,t 2 ] x S as a globally 
hyperbolic manifold with boundary dM = S_|_ U S_. Without reference to the splitting (1), 
M can be characterized as M = J+(S_) n J“(S+) where denotes the causal future and 
past, respectively. Throughout the paper we will assume that n is odd, i.e. the dimension of 
M is even. 



Fig. 1. The Lorentzian manifold X 


1.1. Spinors and the Dirac operator. Let SM —>■ M be the complex spinor bundle on 
M endowed with its invariantly defined indefinite inner product (•,•). Suppose that T> : 
C°°{M; SM) —>■ SM) is the Dirac operator acting on sections of SM. Locally, if 

Co, ei,..., Cn is a Lorentz-orthonormal tangent frame, the Dirac operator is given by 

n 

1=0 

where 7 (X) denotes Clifford multiplication by X, V is the Levi-Civita connection on SM, 
and Ej = {ej,ej) = ±1. Clifford multiplication satisfies 

^{X)^{Y) + -f{Y)^{X) = -2{X,Y) 


and 

(7(X)u,u) = (u,7(X)u) 
for all X,Y £ TpM, u,v £ SpM and p £ M. 

Clifford multiplication with the volume form T = i"'(’^+3)/2 ^(gg) •.. 7 ( 6 ^) satisfies T^ = 
IsM- This induces the eigenspace decomposition SM = S~^M © S~M for the eigenvalues ±1 
into spinors of positive and negative chirality. Since the dimension of M is even, r 7 (X) = 
—7(X)r for all X £ TM. In particular, S^M and S~M have equal rank and the Dirac 
operator anticommutes with T, i.e. it takes the form 


V = 




with respect to the splitting SM = S^M (B S~M. Here D : C°°{M-, S^M) C°°{M; S~M) 
and D : C°°{M; S~M) —>■ C°°{M-, S~^M) are first order differential operators. The subbundles 
S^M are isotropic with respect to the inner product (•, •) of SM. 
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Let V be the past-directed timelike vector field on M with {u, u) = —1 which is perpendic¬ 
ular to all Tit- Then the divergence theorem implies for all u,v € C^{M]SM) 

(2) j {{Vu^v) + {u,T>v)}dY = I {'y{i')u,v)dA — j {'y{v)u,v)dA 

where dV denotes the volume element on M and dA the one on T,f We denote the formal 
adjoint of any linear differential operator 

L : SM) ^ SM) 

with respect to (•, •) by Ll. If u and v are supported in the interior of M, then the boundary 
contribution in (2) vanishes which implies = —T>. 

As is usual in the physics literature, we write /3 = 7 ( 1 ^). Then /3^ = 1 and (•, •) = (/3-, •) 
defines a positive definite inner product on SM provided (•, •) was chosen with the appropriate 
sign. Note that in contrast to (•, •) the inner product (•, •) on SM depends on the choice of 
splitting ( 1 ) which determines v and hence j3. 

1.2. Restriction to hypersurfaces. The restriction of S^M to any slice can be naturally 
identified with the spinor bundle of i.e. = S'S*. On this restriction (•,•) is now 

the natural inner product. Clifford multiplication 7 t(X) on ST,t corresponds to i/ 37 (X) under 
this identification. Since for any X G TSt 

{lt{X)u,v) = (/ 3 i/ 37 (A)n,i;) = {i-i{X)u,v) = -{u,\-i{X)v) 

= -{u,P'yt{X)v) = -{Pu,'yt{X)v) = -{u,'yt{X)v) 

Clifford multiplication on S* is skew-adjoint. Moreover, under this identification /? an- 
ticommutes with 7 t(A) and with the Dirac operator At on Along we have for 

u € C^{M-,S^M) 

(3) Pu = -/3(v, + iAt-^R)u 

where H is the mean curvature of T,t with respect to v. 

If M is a metric product near i.e. the metric takes the form —dt^ + g where 5 is a 
Riemannian metric on S independent of t, then v = —d/dt (because v is past-directed) and 
H = 0. Denoting the Dirac operator on (J^,g) by A and identifying spinors along the t-lines 
by parallel transport, then (3) reduces near to 

Moreover, in the product case (5 commutes with d/dt and it always anticommutes with A, 
hence 

The second order differential operators DD and DD are normally hyperbolic. In the 
product case they have the form -^ + A'^. 
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1.3. The Atiyah-Patodi-Singer spaces. For any subset / C M denote by x/ : M —M 
the characteristic function of I. If I is measurable we denote the corresponding spectral 
projections of At by Pi(t) := and the corresponding subspaces of by 

L|(St;S'M|si) := 'S'M|si)). Note that if I is an interval then Pi{t) is a zero- 

order pseudodifferential operator and it therefore acts on the Sobolev spaces •S'MlsJ 

for any s € M. For s > 1/2, restriction to hypersurfaces is defined and continuous and we can 
therefore introduce the Atiyah-Patodi-Singer spaces as 

//ips(M; 5 +M) := {u^H\M-S+M) \ P[o,oo)(ti)(^|s*J = 0 = P(_oo,01(^2)(n|s*,)} , 

Hips{M-,S-M) := {u € H\M;S-M) \ P(_,o,o]= 0 = P[o,oo)(i2)(«|s*,)} • 

We also denote H%p^{M-SM) := S+M) ® p-M) C H^{M;SM). 

Since /3 anticommutes with At^ it leaves SM) invariant and continuously maps 

S^M) to (M; S'^M). Similarly, we define “anti-Atiyah-Patodi-Singer” spaces 
in such a way that their boundary values are orthogonal to the boundary values of functions 
in FF^pg. This means 

FFfAPs(M;5+M) := {u G H%M-,S+M) \ P(_oo,o)(ti)(«|s*J = 0 = P(o,oo)fe)(tx|s*J} , 
FFfAPs(M;5-M) := {n € H\M;S-M) \ P(o,oo)(ti)(u|E,J = 0 = P(_oo,o)(i 2 )(«|s*J} , 
and FFfAPs(M; SM) := Hl^p^{M- S+M) © FF|aps(M; P-M) C FF^(M; SM). 

2. The Cauchy problem 

We recall that M is foliated by the smooth spacelike Cauchy hypersurfaces where 
t G [^1,^2]- For any s G M the family {FF^(St; ^t 2 ] is a bundle of Hilbert spaces 

over the interval [ti,t 2 ]- As a bundle of topological vector spaces, it can be globally trivialized 
by parallel transport along the t-lines, for example. Continuous sections of this bundle will 
be called spinors of finite s-energy and the space of all such sections will be denoted by 
FE^{M] S^M). Any spinor u of finite s-energy can naturally be considered as a distributional 
spinor on M via 

u[p\= j [ ^^W[(^</^)|sJdt 
J t\ J Sf 

for all test sections ip G {S~^)*M). Note that the volume element of M is AdtdA 

which explains the appearance of the lapse function A as a factor in the above formula. The 
space FE^{M; S~^M) is topologized by the norm 

||'w||ees = max ||u(t)||H'*- 

For s > ^ the Sobolev embedding theorem implies FE^{M] S~^M) C C^{M-, S~^M). 

In order to treat the inhomogeneous Cauchy problem for the Dirac operator we denote the 
space of Fy^-sections of the bundle {F[‘^{T,t', S~^M\Y,t)}te[ti,t 2 ] ^2]; FF®(S,)) and equip 

it with the corresponding L^-norm 

\Ml2^H‘--=[ ||(Au)|sJ|l^«dt. 

Jti 

Now we define 


FE^{M;D) := {u G FE^{M;S+M) \ Du G 
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with the norm 

Here D is applied to u in the distributional sense. The following statement is known as 
well-posedness of the inhomogeneous Cauchy problem for the Dirac equation: 

Theorem 2.1. For any s € M and any t G [^ 1 ,^ 2 ] the mapping 

resi © D : FE^{M- D) ^ W{Et, 5+M|sJ © 

u 1 -^ (rt|st, Du), 

is an isomorphism of Banach spaces. 

Proof. Recall that D = /3 — iAt — fD). Applying Theorem 2.3 in [20, Ch. IV] with 

K = \At + '^H gives that the map is bijective. It is clearly continuous. By the open mapping 
theorem it is an isomorphism. □ 

Remarks 2.2. (a) Since S+Mjst) © is a Hilbert space we con¬ 

clude that FE^{M;D) also carries a Hilbert space topology. Note that FE^{M; S^M) is only 
a Banach space. 

(b) Smooth sections are dense in Ff^CPp, S~^M\Y,t ) and in L'^{[ti,t 2 ]', Since the 

solution of the initial value problem Du = /, njst = uq, is smooth if / and uq are smooth (cf. 
the proof of Theorem 3.5), we get that smooth sections are also dense in EE^{M] D). 

(c) For general s the spaces EE^{M] D) may depend on the choice of the foliation of M into 
Cauchy hyper surfaces. For s = 0 however, the right hand side in Theorem 2.1 is 

5+M|sJ © i2]; T^(S.)) = 5+M|sJ © S-M) 

and hence independent of this choice. Thus FE^{M] D) is the completion of C°^{M; S~^M) 
with respect to a norm which does not depend on the choice of foliation. Hence EE^{M]D) 
is intrinsically defined by M. 

(d) The space EE^(M; D) is the completion of C°°{M; S~^M) with respect to the norm 

II^IIIeo n = max ||u|s, ||i 2 + ||Du |||2 . 

Te[ti,t2] 

By the previous remark the norm 

Ikist 11^2 + ||Du||^2 

for fixed t G [ti,t 2 ] is an equivalent norm on FE^{M; D). This is not surprising because the 
“energy estimate” (2.5) in [20, Ch. IV] ensures existence of a constant C such that 

ll'a|stlli2 < C ■ (l|u|s^||i2 + ||Du||^2 ) 
for any t,r G [ti,t 2 ]. Integration with respect to r yields 

(i2-ii)||M|sJli2 <C-(||u]||2 + (t2-ti)||Du||^2). 

On the other hand, we clearly have [|u ||^2 < C ll'^|s^||^ 2 - Therefore the L^-graph 

norm for D given by ]|u ||^2 + ||Dm ]|^2 is also an equivalent norm on EE^{M; D). 

Both EE^{M-, S~^M) and EE^{M-,D) induce the same relative topology on 
FD^(M;ker(D)) := {u G EE^M'S^M) \ Du = 0}. 

We get well-posedness of the homogeneous Cauchy problem for the Dirac equation: 
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Corollary 2.3. For any t G [^ 1 ,^ 2 ] the restriction mapping 

rest : ker(D)) ^ 5+M|sJ 

is an isomorphism of topological vector spaces. 

For t,t' G [ti,i 2 ] we define the wave evolution operator 

by the commutative diagram 


FE%M;kev{D)) 





Q(t' ,t) 





□ 


2.1. Properties of the evolution operator. It is clear from the definition that for each 
s G M the operator Q{t',t) is an isomorphism of topological vector spaces and that Q{t",t') o 
Q{t',t) = Q{t",t) holds for all t,t',t" G [ti,t 2 ]. In particular, Q{t,t) = Q{t^t') = 


Lemma 2.4. For s = 0 and any t,t' G [ti,t 2 ] the wave evolution operator is unitary, i.e. an 
isometry 

Q{t',t) : L2(St;5+M|s,) ^ 5+M|s„) . 

Proof. It suffices to check ||Q(^^ ,) = ||'^^||L 2 (St) for u in a dense subset of 

Fix s > ^ + 2 and let u G S+MlsJ. By Corollary 2.3 there 

is a unique <1> G FE^{M‘,kei{D)) which restricts to u. Since = 0 equation (3) 

shows Vt4> = -NV^ = N{-iAt + and hence G FE^-\M; S+M). There¬ 
fore <1> is a C^-section of S~^M\Y,t)}te[ti,t 2 ]- By the Sobolev embedding theorem 

C Ci(St;5+M|sJ. Hence 4> G C\M;S+M). 

Thus (2) applies to on x S C M (w.l.o.g. assume t' > t) and gives 

0 = / {(£)$,$) +($,T)$)}dV 
Jm 

= / {fdQ{t',t)u,Q{t',t)u)dA — j (/3m, u)dA 
= \\Q{t ■,t)u\\^2(Y,^,) ~ ll^llL2(Si) • fo 


With respect to the L^-orthogonal splittings 

; 5+MI2,J = {Et, ; (S*, ; J , 

we write Q{t 2 ,ti) as a 2 x 2-matrix, 


(4) 


Q{t2,tl) = 


( Q++(^2Ti) 

I Q_+(t2,tl) 


Hence Q — {t 2 ,ti) = P(-oo,o]it 2 ) o Q{t 2 ,ti)\L 2 
three entries in the matrix. 


—oo,0) 


<3+-(i2Ti)A 

Q—{t2,ti)) 

(St^;5+iW|Et ) ^^*4 similarly for the other 
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Lemma 2.5. The operator Qj^ _ {t 2 ,ti) restricts to an isomorphism 


ker[Q— {t 2 ,ti)] ker[Q++(i 2 , ti)*] 

and Q _ \-{t 2 ,ti) restricts to an isomorphism 


ker[Q++(t 2 ,ti)] ker[(5—(t 2 ,ii)*]- 

Proof. Since Q{t 2 ,ti) is unitary we have Q{t 2 ,ti)*Q{t 2 ,ti) = 1. Spelled out in terms of the 


matrix entries this means 


Q++{t2,tl)*Q++{t2,tl) + Q — \-{t2,tl)*Q-+{t2,ti) = 1, 

Q - {t2,tl)*Q - {t2,tl) + Q+-{t2,tl)*Q+-{t2,ti) = 1, 

Q-+it2,tl)*Q - {t2,tl) + Q++{t2,tl)*Q^ — {t2,tl) = 0, 

Qh — {t2,tl)*Q++{t2,tl) + Q {t2,tl)*Q — |_(t2,tl) = 0. 


(5) 

( 6 ) 

(7) 

( 8 ) 


If n G ker[Q_(^ 2 ,^ 1 )] then by (7) 

Q++{t2,ti)*Q+-{t2,ti)u = -Q-+{t2,ti)*Q — {t2,ti)u = 0. 

Hence Q_|_ {t 2 ,ti)u G ker[Q++(t 2 , ti)*], i-e. _(^ 2 ,^ 1 ) restricts to a map ker[Q_(^ 2 )^ 1 )] —>■ 

ker[(5++(t2,ti)*]. Moreover, by (6), 

U = Q - {t 2 ,ti)*Q - {t 2 ,tl)u + Q+-{t 2 ,tl)*Q+-{t 2 ,ti)u = Q+-{t 2 ,tl)*Q+-{t 2 ,tl)u. 

Similarly, using Q{t 2 ,ti)Q{t 2 ,ti)* = 1, one sees that Q-i_ {t 2 ,ti)* restricts to a 

map ker[Q++(t 2 ,ti)*] ^ ker[(5_ {t 2 ,ti)] and that for u G ker[Q++(t 2 ,ti)*] we have 

Q+-{t 2 ,ti)Q+-{t 2 ,ti)*u = u. Thus Q+-{t 2 ,ti)* : ker[(5++(t2,^i)*] —S' ker[Q—(^ 2 ,^ 1 )] is 
the inverse of Q+-{t 2 ,ti) : ker[Q—(^ 2 ,^ 1 )] ^ ker[( 5 ++(t 2 , 

The proof of the second statement is analogous. □ 

Lemma 2.6. The operators _(^ 2 ,^ 1 ) o P{-oo,o)(ti) Q _ ^{^ 2 ,^ 1 ) ° P{o,oo)iti) 

) continuously to In particular, _(^ 2 ,^ 1 ) cmd 

Q _ ^{t 2 ,ti) are compact as operators from ) —>■ iP'{Ta 2 ’T 

^ 7.2(St2;S’+M|st2), respectively. Similarly, Q_+(t 2 ,ti)* oP(_^ o)(t 2 ) 
and Q_|_(t 2 j^i)* ° P{o,c»D){t 2 ) map 5'+M|st^) continuously to S~^M\-s^^). 

Proof. We show this only for (^ 2 )^ 1 ) and its adjoint. The proof for Q |_(t 2 ,ti) is anal¬ 
ogous. We will need a description of the principal symbol of (5(^2, ^ 2 ) as a Fourier integral 
operator. The operator DD is a normally hyperbolic operator on the globally hyperbolic 
manifold M and, by Theorem A.l, the solution operator 




of the initial value problem 


DDu = 0, ^12^=0, (-Vi.rt)|st =/, 


is therefore a Fourier integral operator of order — | with principal symbol given by (26) and 
canonical relation Ci as described in Theorem A.l. The map D oTto (3 : f ^ u solves the 
Cauchy problem 


Du = 0, tt|st = /, 
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and by the above it is a Fourier integral operator of order — ^ and canonical relation Ci. 
We obtain Q{t 2 ,ti) = rest^ o D oJl^o fi. The restriction operator resjj is a Fourier integral 
operator of order ^ and the canonical relation of resjj is 

^2 = {((y,r?),(x,0) e X T*M I (x,C) € t*M, res,;(x,0 = {v^v)] ■ 

The composition C = C 2 o Ci \s thus described as follows. If (y,T/) is in T*Tit^ then 
there exist two lightlike covectors (y,r/+), € T*M that restrict to {y,r]). For the sake 

of dehniteness we choose fj+ to be future directed and ? 7 _ to be past directed. The orbit 
of {y^V±) under the geodesic flow intersects r*M|sj^ at precisely one point (x±,^±). Let 
(2^±)C±) be the pull back of these co-vectors to T*T,t^. The canonical relation C therefore 
relates {y,r]) to the two points (x+,|+) and (x_,^_). Since M is globally hyperbolic the 
composition C 2 o Ci is proper and transversal (see the discussion in Section 5.1 of [11]). This 
implies that resij o D o Tt^ o /3 is Fourier integral operator of order zero (see for example 
[16, Thm. 25.2.3]) with canonical relation C. In particular, restj o D o Tt^ o /3 is a sum of 
two Fourier integral operators and F-, such that F± is associated to a canonical map 
(y,r/) i-A (x,^±). The principal symbols of these Fourier integral operators can be computed 
using [16, Thm. 25.2.3] or [11, Thm. 4.2.2] and are given by 

(-f±(j^)/? + 7t2(^±)) ° o / 3 , 

where F^^ is the operator of parallel transport from (y, y±) to (x, ^±). Note that the 

projections P(o,oo) (^ 2 ) aud P(-oo,o) (f 1 ) are a pseudodifferential operators of order zero acting on 
L^(St 2 ; ). Their principal symbols Up^ooojlL) ^P(-ooo){h) be computed from 

the principal symbol aAt of the operator At, which is easily found to be (JAtiv) = 

One obtains 

^P(0,oo)it2)iy^v) = \{^- /3|l^irSt2(^)) > 

^P(-^fi){ti){y^y) = \ ( 1 + /3||^rSti(7)) • 

Therefore the operator _(^ 2 ,^ 1 ) o ^’(-oo,o)(fi) is a sum of two Fourier integral operators 

associated to the canonical maps (y, y) i-A (x,^±)with principal symbols 

g±(y,y) = ±^||^±ir^ (||^±|| - ht 2 {i±)) (-f±(i^)^ + 7t2(C±)) 

°^(^,C±),0,r;±) ° Il^ir^/3(ll7ll +/3 Ti(7))- 

Since is lightlike we have ^±(z^) = ±||^±||. Thus 

(lie+ll - ht2{i+)){-U{^)li + 7t2(e+)) = - (lie+ll - ht2{i+)) (IIC+II + /37i2(e+)) /? = 0, 

and consequently q+{y,r]) = 0. We have used here that /3 anticommutes with 7 * 2 (C+)- Clif¬ 
ford multiplication by is given by (t||'^±||/3 + 7 t 2 (^±)) and because parallel transport is 
compatible with Clifford multiplication we have 

r(a;,i±),(y,i)±) ° + = (Til^eill/s + T2(?±)) 

Thus 

q-{y,v) = -^ll?-ir^(^ll?-ll +7t2(^-)) o \\r]\\-^ (||y|| 
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= o/ 3(||7?|| +/37ti(?7)) (M - Pjtiiv)) ^ 

= 0 . 

We have shown that the principal symbols of the zero order operators vanish and therefore 

_(^2) ii)o-P(-oo, 0 )(^i) is the sum of two Fourier integral operators of order —1, each of them 

associated to a canonical map. The statement of the theorem now follows from the mapping 
properties of Fourier integral operators whose canonical relation is a canonical graph (see for 
example [16, Cor. 25.3.2]). □ 

Corollary 2.7. The kernel of Q _ and Q++{t 2 ,ti) respectively, consists of smooth 

sections. 

Proof. Equation (6) implies that every element in the kernel of Q _ {t 2 ,ti) is in the (+1)- 

eigenspace of _ {t 2 ,ti)*Q^ _ {t 2 ,ti). Lemma 2.6 then shows that this eigenspace consists of 

functions with Sobolev regularity of any degree. The same argument applies to the kernel of 
Q++it2,ti). □ 

3. The Dirac operator with Atiyah-Patodi-Singer boundary conditions 

In this section we relate the Fredholm property of Q _ {t 2 ,ti) to that of the Dirac oper¬ 

ator itself under Atiyah-Patodi-Singer boundary conditions. The discussion is based on the 
following auxiliary lemma. For a proof see e.g. [4, Prop. A.lj. 

Lemma 3.1. Let H be a Hilbert space, let E and F be Banach spaces and let L : H ^ E 
and P : H ^ F be bounded linear maps. We assume that P : H ^ F is onto. 

Then L|j^gr[P] • ker[P] E is Fredholm of index k if and only if L (B P ■ H ^ E (B F is 
Fredholm of index k. □ 

Theorem 3.2. The operator 

(^[0,oo)(ii) oresq) © (P( 

— CXD,0 ](t 2 ) o resq) ® D : 

FE\M- D) ^ © L\M- p-M) 

is Fredholm and its index satisfies 

ind[(P[o,oo)(^i) o resq) © (P(-oo, 0 ](^ 2 ) o resq) © D] = ind[(5—(^ 2 ,^ 1 )] • 

Proof. Theorem 2.1 with s = 0 shows that FE^{M-, D) carries the topology of a Hilbert space 
and that 

D : FE\M;D) ^ L‘^{[ti,t2], L\E,)) = L\M;S-M) 
is onto. We can therefore apply Lemma 3.1 with 
H = FE^{M;D), 

E = ^[0,oo)(Sq; 

F = L^{M-S-M), 

L = {P[o,oo){ti) o resq) © (P(_oo,o](^ 2 ) o resq) , and 
P = D. 

It remains to check that the operator Lq := L|ker(D) • FE^{M;kei{D)) E is Fredholm with 
the same index as Q _ {t 2 ,ti). As to the kernel we have 

ker[Lo] = {ue PE0(M; ker(D)) | P(-oo, 0 ] (t2)(R|Eq) = 0 } 
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= ker[Q —(t 2 ,ii)] • 

The image is given by 

im(Lo) = {{ui,U 2 ) I U 2 = Q-+{t 2 ,ti)ui +Q —for some w € 

= {{ui,Q-+{t 2 ,ti)ui + Q — {t 2 ,ti)w) I Ui G W G ^^_oo, 0 )(^il)} • 

This shows in particular that the image is closed. Namely, let ui^i G ) 

and Wi G 5+M|stj) such that (upi, Q_+(t2, ^ 1 )^ 1 ,* + Q— {t 2 ,ti)wi) (mi,M 2 ) 

in E. Then Q_+(t 2 , Q-+{t 2 ,ti)ui and hence Q — {t 2 -,ti)wi U 2 - Q-+{t 2 .,ti)ui. 

Since Q _(t 2 Ti) is Fredholm its image is closed and thus U 2 — Q _ \-{t 2 ,ti)ui = Q _ {t 2 ,ti)w 

for some ic G )• Therefore ( 111 , 112 ) = {ui,Q-^{t 2 ,ti)ui + Q —(t 2 ,ii)ii’) 

lies in im(Lo). 

Moreover, 

im(Lo)-^ = {{vi,V2) G E \ (m, 111 ) 2,2 + ( 112 , (5-+(t2, ii)iii + Q — {t2,ti)w)L2 = 0 

for all m G 11 G 

= {{vi,V2) ^E \ (m + Q-+{t2,ti)*V2,Ui)L2 + {Q - {t2,ti)*V2,w)L2 = 0 

for all 111 € w G F(_oo,o)(^ti)} 

= {(i’i,i’ 2 ) G Fi I ui + Q-+{t 2 ,ti)*V 2 = 0 and Q — {t 2 ,ti)*V 2 = 0} 

= {{-Q-+{t 2 ,ti)*V 2 ,V 2 ) I V 2 G ker[Q—(t 2 ,ti)*]} 

^ker[Q__(t 2 ,ti)*]. □ 


We now consider finite energy spinors which satisfy the Atiyah-Patodi-Singer boundary 
conditions and define 

FEXpg{M-, D) := {u G EE''{M]D) \ F’[o,oo)(ii)(^^|st J = 0 = F’(-oo,o] (i 2 )(ii|Et 2 )} • 

The space EE^pg{M] D) is a closed subspace of FE^{M] D) and carries the relative topology. 
Theorem 3.3. The operator 

Daps ■■= : FElp^{M-D) -a L\M;S-M) 

is Fredholm and its index satisfies 

indpAPs] = md[Q—(^ 2 ,^ 1 )] • 

Proof. This time we apply Lemma 3.1 with 
H = FE^{M;D), 

E = L‘^{M;S-M), 

F = (Si,; 5+M|s,,) © i S+M\,,,^), 

L = D, and 

P = {P%oo){ti) o resi,) © (P(-oo,0](^2) o restj) • 

Then we get that Haps is Fredholm with ind[HAPs] = ind[H © P], Theorem 3.2 concludes 
the proof. □ 
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Looking at finite energy spinors satisfying anti-Atiyah-Patodi-Singer boundary conditions 
D) ■= ^ FE^{M-,D) I P(_oo,o)(ii)(«|stJ = 0 = P{o,oo)(f2)(R|st2)} 

the same reasoning shows 
Theorem 3.4. The operators 

{Pi —oo,0 )(ti) O restj © (P(0,oo)(i2) o resjj) © D : 

FE\M- D) ^ S +© L\M-, S-M) 

and 

D.APS ■■= ^ ^^aAPs(M; D) -A L\M; S-M) 

are Fredholm and their indices satisfy 

ind[(P[o,oo)(ii) orestj © (^(-oo, 0 ](^ 2 ) orestj) © L>] = ind[DaAPs] = ind[(3++(t2,ti)] • □ 

In the Riemannian setting Dirac operators are elliptic and hence the kernels of Dirac oper¬ 
ators consist of smooth spinors, independently of any boundary conditions. In our Lorentzian 
situation this is no longer true. Solutions to the Dirac equation can have low regularity. 
Remarkably, under APS or aAPS boundary conditions solutions are again smooth just as if 
we had elliptic regularity theory at our disposal. 

Theorem 3.5. The kernels of the operators Daps ■ FE^pg{M ■, S~^M) S~M) and 

.DaAPS : M) —>■ L‘^{M; S~M) consist of smooth spinors. The indices satisfy 

ind[DAPs] = -ind[DaAPs] = dimker[D|c“pg] - dimkerplc^pj. 

Proof. It was shown in the proof of Theorem 3.2 that 

ker[L)APs] = ker[Lo] = {u G FE^{M-,D) \ Du = 0 and G ker[(5_(^ 2 ,^ 1 )]} ■ 

Since u\s^^ is smooth by Corollary 2.7, the solution u of the Dirac equation lies in FE^(M] D) 
for any s G M by Theorem 2.1. The Dirac equation 

, , du , n 

(9) — = lAtu + -Hu 

implies that u is a C^-section of Differentiating (9) and iterating the argument 

shows that u is a C^-section of for any /c G N and any s G M. By the Sobolev 

embedding theorem u is smooth. The argument for the kernel of DaAPS is the same using 
Q++{t 2 ,ti) instead of Q — {t 2 ,ti). 

Therefore we have ker[DAPs] = ker[D|c'|°pg] and, by Lemma 2.5, 

dim coker [Daps] = dimker[Q—(^ 2 ,^ 1 )*] = dimker[(3++(t2Ti)] = dimker[D|c^pg] . 

The argument for ind[DaAPs] is analogous. □ 

4. The geometric index formula 

In this section we give a geometric expression for the index considered in the previous 
section. We use the splitting M = [fi,f 2 ] x S and the Lorentzian metric —N'^dt^ + gt to 
introduce the auxiliary “Wick rotated” Riemannian metric g := pN'^dt^ + gt- The formula 
analogous to (3) for the Riemannian Dirac operator T) reads 

(10) Dn =-/3 (-V, + At + 


u. 
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Since the induced metrics gt on the slices have not changed, the operators in (3) and 
in (10) coincide. Now we deform g near the boundary dM to a metric g so that g = dt^ + gt^ 
near {ti} x S and g = dt'^ + gt^ near {^ 2 } x S. Then the corresponding Dirac operator 
i) = + Bt) : C°°(M; 5+M) ^ S-M) simplifies to 




near {tj} x S where j = 1,2. 

For any subset / C M denote by x/ : M —)■ M the characteristic function of I. Denote the 
corresponding spuectral projections of At by Pi{t) = xii^t) and the corresponding subspaces 
of by 5+M|sJ := S+M\j:,)). 

The operator t) : S~^M) —>■ S~M) is Fredholm where 

(11) i7lps(M;5+M) = {uG H\M-,S+M) \ P[o,oo)(ii)(«|s,J = 0 = P(-oo,o](i 2 )(n|s,J} • 

By the Atiyah-Patodi-Singer index theorem [1, Thm. 3.10] the index is given by 


( 12 ) 


ind[Z)] 



h(Afi) + g{At^) + h{—At2) + r/(—Af^) 
2 

h{At ^) + h{At^) + g{At^) - g{At^) 

2 


where A(V) is the A-form computed from the curvature of the Levi-Civita connection V of 
g and where h{A) denotes the dimension of the kernel of A and g{A) its ry-invariant. 

Next we express the A-integral in terms of the Levi-Civita connection V of the original 
Lorentzian metric g rather than that of g. Let a be the AdQL(n+i,R)-invariant polynomial 
giving rise to the A-form, 

A(V) = a(D^,...,D^). 

Here denotes the curvature 2-form matrix of V with respect to any local frame of TM. 
Denote the connection 1-form matrix of V by io^. We define the trangression form by 

(13) TA(V, V) / a(w^ 

2 Jo 

It is well known that 

A(V) -A(V) = (iTA(V,V). 

The Stokes theorem yields 

(14) / A(V)=/ A(V)+/ TA(V,V). 

Jm Jm JdM 

Equation (13) implies easily that the pullback of TA(V, V) to dM is a polynomial expression 
in the curvature tensor of V, the second fundamental form of the boundary with respect to 
V, its first derivatives, and the first and second derivatives of the lapse function N. Thus 
TA(V,V) is determined by the geometry of {M,g) and does not depend on the particular 
choice of auxiliary metric g. Henceforth we write TA(( 7 ) instead of TA(V, V). Inserting (14) 
into (12) yields 


ind[D] = [ A(V) + [ 

Jm JdM 


h{At^) + h{At2) + g{At-f) — g{At^) 


(15) 


TA(g) 


2 
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4.1. Spectral flow. We put := ^{h{A) + ri{A)). Then the spectral flow of the operator 
family {At)tx<t<t 2 is defined as 

sf(Ae[ii,t2]) = 3 {t2) 

where £,{Af) = c{t)+ jif) is the splitting into a continuous function c{f) and an integer-valued 
function j{f) with j{ti) = 0, cf. [2, Sec. 7]. We rewrite (15) as 

ind[.D] = [ A(V) + [ TA{g) + - C{At,) - h{At,) 

Jm JdM 

= [ A(V)-h / TA{g) + c{t 2 ) + j{t 2 ) - c{ti) - h{At^) 

Jm JdM 

We rearrange the terms as 

(16) md[D] - j{t2) + hiAt^) = [ A(V) + / TA(c/)-h c(t2) - c(ti). 

Jm JdM 

All terms on the left hand side of (16) are integer valued. On the other hand, if we consider 
the right hand side 


/ 

J\t 


[tl,l2]xS 


A(V) + 



TA(<7) 



TA( 5 ) + c{t 2 ) 


c{tl) 


as a function of t 2 while keeping ti fixed, then all terms are continuous in ^ 2 - Hence both 
sides in (16) are constant. In the limit t 2 —>■ the RHS tends to 0. Therefore 


(17) ind[T»] - + ^(^* 2 ) = / A(V) [ TA{g) + c{t 2 ) - c(ti) = 0 . 

Jm JdM 

We conclude 


(18) 


sf(^te[ti,t2]) 


ind[Z)] + h{At^) 

[ A(V) + f 
Jm JdM 

[ A(V) + f 
Jm JdM 


) - h{At2 ) + vjAti ) - g{At ^) 
TA{g)+C{At2)-C{At,). 


4.2. Fredholm index of the wave propagator. Next we relate ind[(5 (^ 21 ^ 1 )] to the 

spectral flow sf(Aig[i^ ^j])- We fix ti and consider ind[Q (^Ti)] as a function of t. For 

t = ti, the operator Q _(ti,ti) is the embedding 

A?_oo,0)) 

and hence has 


(19) 


ind[(5—(ti,ti)] = -h{Atf). 


We use the characterization of spectral flow described in [19]. We choose a partition ti = 
tq < Ti < ■ ■ ■ < tn = t 2 and numbers Oj > 0 such that ±aj ^ spec(At) for all t € [Tj-i,Tj]. 
Then 


N 

(20) sf(A,g[i^,,,]) = ^ (dimLfo^„^.](S.^.;5+M|s.p -dimLfo_,^,](S.^_,;5^ 
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Since aj ^ spec(^t) the family of projectors -P(_oo,aj] (^) is a continuous section of bounded 
operators on the Hilbert space bundle S~^M\Y,t)}t over [Tj-i,Tj]. Thus 

Pi-oo,a,]it) o Qit,h) : ^ 

is a continuous family of Fredholm operators for t E [Tj-i,Tj]. In particular, 

(21) ind[P(_oo_„.](r,) o Q{Tj,ti)] = ind[P(_oo,a,.](H-i) ° Qi^j-iPi)]- 

If we consider both operators P(-oo,a-]{'t) ° Qi^Pi)) aiid P{-oofl]{'t) ° QitPi)) as operators 
P‘{-oo,0)(P‘ti'^ ^(_oo,a,](^d‘S'’^^|st), then they differ by P(o,a,■](i)oQ(i,H))- Since 

the spectrum of Ar is discrete, this difference operator is of finite rank and hence compact. 
Therefore 

(22) ind[P(_oo,a,](i) « QitPi))] = ind[P(_oo,o](i) o QitPi))] 

where both operators are considered as operators from ) to 

Now Q _is the same as T’(_oo,o](0 ° Qi'tPi))i except that it 

is considered as an operator to gj (S^; 5+M|st). Hence 

(23) ind[Q__(t,ti)] = ind[P(_oo,o](0 « <3(^> ^i))] + dimL(o^„.](St; 5+M|st) 
for t E [Tj-i,Tj]. We conclude 

N 

ind[Q— it 2 ,ti)] = ^(ind[Q— (rjAi)] - ind[Q—(rj_i, H)]) + ind[Q—(ti,ti)] 
i=i 

N 

(^ind[P(_oo,o](H) ° Qi^jPi)] + dim^?o,a,](^rH ‘S+M|s,p 
i=i 

- ind[P(_oo,o](H-i) ° 

+ ind[Q— 

N 

= Y (ind[P(-oo,0](H) ° Qi'TjPi)] + dimLfo ,,.](S^.; S’+M|s,.) 
i=i 

- ind[P(_oo,o](H-i) o Q{Tj-i,ti)] - dimLjj ,^^](S^^_^; S’+M|s,,_J 

- h{Ar^) + h{Ar^_^)J + ind[Q—(Hpi)] 

N 

= Y (™d[P{-oo, 0 ](p) o Q{rj,ti)] - ind[P(_oc^ 0 ](H-i) o <3(p-iPi)]) 

i=i 

+ sf(^te[ti,i 2 ]) “ ^{At^) + h{Atj^) + ind[(5—(HPi)] 

N 

= Y (ind[P(-oo,0](H) ° Qi^jAi)] - ind[P(_oo,o](p-i) ° 

i=i 

+ sf(^tg[ti,t2]) — h{At^) 
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N 

= ^ (^ind[P(_oo,a,]('ri) o Qirj,ti)] - ind[P(_oo,a,]('ri-i) o QiTj-i,ti)]^ 
i=i 

+ — h{At^) 

[ A(V) + [ TA{g) - ^(^^ 1 ) + ) - r]iAt2 ) 

Jm JdM 2 

We have proved: 

Theorem 4 . 1 . The index of Q —(^2,^1) : -^(-oo,o)(^ii’^ ^\T.t2) 

is given by 

ind[Q__(t 2 ,t.)]= / A(V)+ / + ^ ^ 

Jm JdM 2 

5. An example 

We describe an example with nontrivial index in any dimension divisible by 4. Let A: G N. 
The manifold will be M = [ti, t 2 ] x equipped with a metric of the form —dt^ + gt where 

gt is a 1-parameter family of metrics on the sphere arising as follows: 

Consider the Hopf hbration ^ If carries its canonical metric of 

constant sectional curvature 1 and the Fubini-Study metric with sectional curvature 

between 1 and 4, then the Hopf hbration is a Riemannian submersion. The hbers are great 
circles. We now rescale the standard metric of by the factor t > 0 along the hbers 

and keep the metric unchanged on the orthogonal complement to the hbers. This yields a 
1-parameter family of metrics gt on known as Berger metrics. 

All Berger metrics are homogeneous under the unitary group U(2/c). Representation the¬ 
oretic methods have been used to compute the full spectrum of the Dirac operator At on 
(54^-1, grj) in [14, Prop. 3.2] for the case k = 1 and in [3, Thm. 3.1] for general k. As a 
corollary one obtains a Dirac eigenvalue X{t) = (—l)^“4^i _ 2 k) of multiplicity (^). This 
eigenvalue vanishes for t = Ak while all other Dirac eigenvalues on {S^^~^,g 4 k) are nonzero. 
Hence, if we choose ti slightly smaller than 4A; and t2 slightly bigger, then the spectral how 
is sf{At^^tiM') = We conclude 

ind[DAPs] = ind[Q__(t2,ti)] = sf(Ag[q,t2]) = ■ 

6. A parametrix for the APS Dirac operator in case the metric near the 

BOUNDARY IS OF PRODUCT TYPE 

Since the Dirac operator subject to APS boundary conditions is a Fredholm operator its 
equivalence class in the Calkin algebra is invertible. In the special case when the Lorentzian 
metric is of product type near the boundary components it is in fact possible to construct 
an inverse modulo smoothing operators. 
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Fig. 2. The manifold M with product structure near the boundary 
Such a parametrix has the same microlocal properties as the Feynman propagator and we 
will refer to it as a Feynman parametrix. The Feynman propagator is a fundamental solution 
of the Dirac operator that was introduced by Richard Feynman to describe the physics of 
the electron in Quantum Electrodynamics. Feynman parametrices appear in quantum field 
theory on curved spacetimes and were constructed as one type of distinguished parametrix in 
the seminal paper of Duistermaat and Hormander [12]. 

6.1. Feynman parametrices. Let L be a normally hyperbolic operator on X acting on the 
sections of a vector bundle E, see Appendix A for details. Let G be a left parametrix for L, 
i.e. G is a continuous map G : G^(A; E) C°°{E) and 

GL = idc^ + R, 

where R is an operator with smooth integral kernel. As usual, define WF^(G) as the set 

WF'(G) = G t*{X X X) | G WF(G)}, 

where T*{X x X) = T*{X x X)\0, and WF(G) is the wavefront set of the distributional 
integral kernel of G. Let A* be the diagonal in T*{X x X) and <I>i the geodesic flow on T*X. 
We put 

Kp := A* U {(x, Q x'€ T*{X x X) \ ^ is lightlike and 31 > 0 : ^t{x', ^') = (x, ^)}. 
Then we say that G is a Feynman parametrix if 

WF'(G) C Af- 

Feynman parametrices are known to exist and are unique up to smoothing operators on 
any globally hyperbolic spacetime (see [12, Thm. 6.5.3] taking into account that globally hy¬ 
perbolic spacetimes are pseudoconvex with respect to normally hyperbolic operators). Every 
Feynman parametrix is also a right parametrix, i.e. 

LG = idc^ -|- R, 

and extends to a continuous map G : Hq{X;E) —>■ HfJ^^{X-, E) for all s G M. 

Analogously, a left parametrix S for a Dirac type Dirac operator T> called a Feynman 
parametrix if WFQ5) C Ap- By the same reasoning as in [12, Thm. 6.5.3] such parametrices 
are unique up to smoothing operators. Since L = is normally hyperbolic there exists a 
Feynman parametrix G for L. The operator is a direct sum of DD and DD. The Feynman 
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parametrix G can therefore also be assumed to be the direct sum of Feynman parametrices 
for DD and DD, respectively. 

Lemma 6.1. The right parametrix S := VG for V is a Feynman parametrix. It maps 
Hq{X] S^X) continuously to H^^^{X-, X) for all s € M. 

Proof. The wavefront set relation is preserved by the composition with V from the left. The 
only thing to show here is that 5 is a left parametrix. This follows from a variation of the 
argument for the uniqueness of distinguished parametrices in [12]. Namely, let S' be the left 
parametrix GV, and define R := S — S' . We will show that R is smoothing, i.e. WF'(R) = 0. 

Clearly, WF'(i?) C Ap, so let G Ap. This implies that either (x,^) = 

or there exists a lightlike geodesic connecting (x,^) and (x',^'). Denote this geodesic arc by 
K G X and choose (p G G^(X) such that (p(x) = 1 for all x in an open neighborhood of K. 
Then, since the first order operator pV — Vp is supported away from K, the point (x, x', ^') 

is not in the wavefront set of (pT> — Rep. Since 

S{p'D — T>p)S' = Sp — pS' 

the point (x,^;x',^') is not in the wavefront set of Sp — pS' . Since this is true for any 
(x, x', ^') G Ap the operator S — S' has smooth kernel. □ 


6.2. A Feynman parametrix in the prodnet case. It is instructive to see what a Feynman 
parametrix for the normally hyperbolic operator L = T>^ looks like in the case A = M x S 
equipped with a product metric. In this case L has the form 

The operator |A| ^ as given by spectral calculus has domain ker[A]“’- so that the operator 
|A|“^(1 — Pq) is well defined and bounded on LfiYi;SM). Here Pq = Pjo} is the orthogonal 
projection onto the kernel of A. Then the family of operators 

ait) ■■= (e"''^'*X(o,oo)(i) + e‘l^l*X(_oo,o)W) l^r^(l - -Po) + i X(o,oo)(i)^o 

dehnes a Feynman parametrix via 

{Gu)it,-)= [ git - s)uis,-)ds. 

Jm. 

For any open interval I the restriction of G to / x S then defines a Feynman parametrix on 
the globally hyperbolic manifold / x S. 

Assume s > 1/2. For M C A of the form [ti,t2] x S we have for any u G PPiM] SM) the 
following formulae, easily obtained by integration by parts: 


(24) 

and 


d 


dt 


d 


—— iAG — + iA]u = u + Ru, 


dt 


d 


dt 


d 


G —±iA - —±iAGU = Qrt 


dt 


Here the map R : H^iM; SM) —)■ H'^ SM) is defined as 
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The operator Q : H^{M; SM) —>■ SM) is given by 


{Qu){t,-) :=--|yi|-i(i-Po 


,-i|A|(i-ii) 




_ pAA\{t-t2) 


{u\ 


+ {t- ti)Po{u\^t ). 


We define the Feynman propagator as S := T>G. One computes directly that 
{Su){t,-) = 13 - i Gu^ (t, •) = (3 J^hit - s)u{s,-)ds , 

where 

Ht) = ^(0,oo)We"‘'^'*X(0,oo)(i) - -P(-oo,0)We''^'*X(-oo,0)(0 + X(o,oo)(i)^o • 

Thus the Feynman propagator propagates elements in the positive spectral subspace of A 
forward in time and elements in the negative spectral subspace of A backwards in time. 

Theorem 6.2. Let M = [ti,t 2 ] x S 6e equipped with a product metric —dt^ + g and let 
S = T>G he the Feynman propagator. Let s > 1/2. Then we have 

•SF\h’>{M;SM) = id//s(M;SM) + R- 

Moreover, S maps H^{M] S~M) to sections satisfying AFS boundary conditions at ti, i.e. if 
u G H^{M; S~M) then P[o^oo)(^i)((‘5ri)|sti) = 0- Moreover, 


— id 


HXpJM-,S+M)- 


Proof. That FS\h‘>{M;SM) = idH»(M;SM) can be checked by direct computation and 
^'^\h‘’{m-,sm) = id/fs(M;S'M)+ -R follows immediately from (24). Since h{t) = —T’(-oo,o)(^)c^^"^^* 
for t < 0 we obtain Pjo,oo)(^i)((‘5tt)|Et^) =0. If ti € H^{M; S~M) then Su G 
and this means it satishes APS-boundary conditions. The formula = 

idH],^pg{M-,s+M) follows from the fact that the operator R vanishes on S~^M). □ 

6.3. Gluing of Feynman parametrices. Let {X, g) be a globally hyperbolic manifold dif- 
feomorphic to (—T, T) x S such that {t} x S is a Cauchy hypersurface for each t G (—T, T). 
Let 0 < (5 < T. The sets Ui = (-T, 6 ) x S, C/a = {S, T) x S and C/ = C/i n C/a = {- 6 , 6 ) x S. 
Let L be a Dirac type operator or a normally hyperbolic operator. Now suppose that Gi is 
a Feynman parametrix on C/i and Ga a Feynman parametrix on C/a. We choose the following 
gluing functions , 0a, </>i, </>2 G G°°{X) satisfying the following properties 

• supp(0i) C supp((/>i) C C/i, 

• supp( 0 a) C supp((/ 2 ) C C/a, 

• 01 + 02 = 1, 

• = 1 for all X in some open neighborhood of supp(0i), 

• <f 2 {x) = 1 for all X in some open neighborhood of supp(02). 

This implies that the differential operator [(fi,L\ = (fiL — Lcfi is compactly supported in 
X, and for all x G supp([L, ijii]) we have 0i(x) = 0, 62 {x) = 1, 4 > 2 {x) = 1. The analogous 
statement holds for [L,(j)i]. 

Lemma 6.3. The operator 

G := 0iGiijii + 02G2</>2 ~ GiGi[(pii L]G 2 (t >2 ~ 02G2[</>2, L]Gi(pi. 
is a Feynman parametrix. 
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Proof. It is easy to see that G is well defined. We first check that WF'(G) C Ap. Since Gi 
and G 2 are Feynman parametrices both WF'(0iGi(?!)i) and WF'{ 62 G 24 > 2 ) are subsets of A^. 
Moreover, by the mapping properties of the canonical relation Ap, we have 

WF'(02G2[02, C €Ap\^€Vp,x€ U2,x' € U^}, 

where V+ C T*M is the closed forward light cone, i.e. the set of future-directed causal covec¬ 
tors. We have used here that the past of supp([(^ 2 , A]) has empty intersection with supp( 02 ). 
Similarly, using that the future of supp([(/)i, L]) has empty intersection with supp(0i), one 
obtains 

WF'(0iGi[(/)!,L]G 2 </) 2 ) C {(x,e;x',e') € I € y+,x G Ui,x' G C/ 2 }. 

This shows that WF^(G) C Ap. Now we check that G is a left parametrix. Let us introduce 
the relation A^B A — B is a smoothing operator. Note that 

0 iGi(/)iL~ 0 i 0 i + 0 iGi[,^i,L] 

and 

(^iGi[4>i, L]G24>2L ~ 0iGi[0i, L]G2[<^2, A] -|- 9iGi[(j)i, L]. 

Collecting these two terms we obtain 

OiGicpiL — 6iGi[4>i, L\G24>2L ~ 0i4>i + Ai, 

where Ri = —OiGiQi, and Qi = [(fi, L]G 2 \(t> 2 ^ A]. Since WF'(G 2 ) C Ap we have 

WF'(Qi) C G AiT’ I X G supp([(/>i. A]),x' G supp([ 02 , A]),^ G V+}. 

Again, the future of supp([(/)!, A]) does not intersect supp(0i) and hence we obtain WF'(i?i) = 
0. Thus, Ri is a smoothing operator. In the same way one shows that 

92G2(t>2L = + A 2 , 

where R 2 is a smoothing operator. All together we obtain 

GL = 0i<^i -|- 0202 + Ri + R 2 = idcg^ + Ri 

where R = Ri + R 2 is a smoothing operator, and hence G is a left parametrix. □ 

6.4. A distinguished Feynman parametrix if AP has product structure near the 
boundary. In this subsection we will be assuming that M is isometric to a Lorentzian 
cylinder near the boundaries and . This means near these boundaries the metric is of 
the form —dt^ + g where (7 is a Riemannian metric on S independent of t. More precisely, 
we assume product structure on -|- e) x S and {t 2 ,t 2 — e) x S. Thus, we can think of 

M as the union of the two product manifolds -|- e] x Ei, [t 2 T 2 — e] x E 2 and the open 

globally hyperbolic manifold X = (ti, t 2 ) x S. Starting with any Feynman parametrix for the 
operator T> on A we can use the Feynman parametrices on -|- e] x Si, [^ 25/2 — e] x S 2 

to obtain a Feynman parametrix S by gluing. The following theorem is immediately deduced 
from the support properties of the gluing functions and Theorem 6.2. 

Theorem 6.4. Assume that the metric on M has product structure near Si and S 2 . Let 
s > 1/2. Then the operator S continuously maps H^{M; S~M)) to S~^M) and we 

have 

L)S\hs(^M-,S-M) = + Qi 

AA)|i^^ps(M;S+M) = id_H-Spg(M;S+M) + A- 
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where Q and TZ have integral kernels that are smooth up to dM. Moreover, the integral kernel 
of Q is compactly supported in {M\dM) x M, and the integral kernel of TZ satisfies APS 
boundary conditions in the first variable. Therefore, for Q and TZ have the following mapping 
properties: 

Q : H^{M- S-M) C^{M\dM-, S-M), 

TZ : H^{M; 5+M) ^ C^psi^; 5+M), 

for any r € M. □ 

The above can be understood as a refined version of the Fredholm property of the operator 
D subject to APS-boundary conditions and it gives an alternative approach to the index 
problem via parametrices. An example is the following refinement of Lemma 2.6 in the case 
of product structure of the metric near the boundary. 

Theorem 6.5. Assume that the metric on M has product structure near and Then 
the operators _ {t 2 ,ti)°T^{-oo, 0 )i'ti) andQ _h(^ 2 ) ii)o.P(o,oo)(^i) have smooth integral kernels. 

Proof. In the proof we will assume that the Feynman parametrix was constructed by gluing 
on a slightly larger spacetime M = x S, where M = [ti,t 2 ] x S and we assume the 

metric to be of product type on Z = [t 2 , x E. We will also assume that the gluing functions 
used in the construction have derivatives supported away from Z. We will need a smooth 
function 0 on M which is equal to one on M and which vanishes near t = t^. We choose 9 
such that on Z it depends on the t-variable only (i.e. such that 9 is constant on E*). Now let 
u G ). Then there exists a unique solution <I> G FE^{M,S~^M) 

of the Cauchy problem = 0, rest^(<h) = 0. The spinor D(0d>) = I3{dt9)^ is compactly 
supported in Z. Since $ satisfies APS-boundary conditions on M we have 

S{P{dt9)^) = 9^ + K^, 

where K is an operator with smooth integral kernel. This shows that Q{t 2 ,ti) equals to 
(restj o S){l3{dt9)^) + (restj o K)^. On Z we decompose ^{t, •) = •) -|- $-(t, •) + ‘ho 

according to the spectral decomposition into positive, negative, and zero spectral subspaces of 

A. Hence _ {t 2 ,ti){u) = (restj o5)(/3(9t0)(<h_-|-‘ho))-|-(restj o A")*!*. Since (3 anticommutes 

with A and dt9 depends only on t, the section (/3(9t^)‘h_)(t, ■) is in the positive spectral 
subspace of A for each t. By the propagation properties of <S and the support properties 

of the gluing functions we obtain (resjj o S){/3{dt9)^-) = 0 and therefore _ {t 2 ,ti)u = 

(restj o S){l3{dt9)^o) + (resjj o A')<h. The map 

u (-> (restj o S){j3{dt9)^Q) 

has smooth integral kernel because the kernel of A is finite dimensional and consists of smooth 
sections. A similar argument works for Q _ ^{t 2 ,ti) o P’(o,oo)(H)- D 

7. Concluding remarks 

To be definite we have restricted our presentation to the Dirac operator acting on spinor 
fields. Minor adaptions to the proofs yield the corresponding results for twisted Dirac op¬ 
erators acting on sections of SM E where A —>■ M is a Hermitian vector bundle with 
compatible connection V^. One only has to replace A(V) by A(V) A ch(V®) where ch is 
the Chern character, TA( 5 () by TA(g() A ch(V^), and the boundary Dirac operators by the 
corresponding twisted Dirac operators on the boundary. This covers the physically relevant 
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cases of hermitian bundles induced by compact gauge groups. For example twisting with a 
complex line bundle allows one to include electromagnetic potentials. 

A word of caution here: If the twist bundle E is a natural bundle induced by the geometry 
of M, e.g. E = TM, then its induced metric is not in general definite. If this happens 
the above remark does not apply. Indeed, it is not clear if the twisted Dirac operator is 
a symmetric hyperbolic system in this case and the whole analysis may require a careful 
reexamination. 

As in elliptic index theory there are many directions in which a generalization should 
be possible. It is conceivable that one could work out an equivariant version, a family ver¬ 
sion, a spatially L^-version or other versions of our index theorem for noncompact Cauchy 
hyper surfaces. The latter may be very useful on physical grounds. 

For the Dirac operator on a compact Riemannian spin manifold with boundary one cannot 
use anti-Atiyah-Patodi-Singer boundary conditions instead of APS-boundary conditions. This 
would not give rise to a Fredholm operator. In the Lorentzian setting however anti-APS 
conditions are equally good. This difference in qualitative behavior can be understood as 
follows: Solving the Dirac equation on a Riemannian product manifold under APS-boundary 
conditions is like solving a heat equation which is possible forward in time but not backwards 
in time. In the Lorentzian setting it corresponds to solving a wave equation which can be 
done both forward and backwards in time. 

As a final remark we would like to note that APS boundary conditions make sense also for 
distributional sections of lower regularity. For example, for the restriction of a distributional 
spinor u to to be well defined it is sufficient to assume that its Sobolev wavefront set 
WFs(u) is contained in the light cone for some s > 1/2. This is automatically satisfied if u 
solves the Dirac equation. Distributional solutions that satisfy APS boundary conditions can 
be shown to be smooth, and therefore allowing lower regularity results in the same theory. 

Appendix A. Solution of the Cauchy problem for normally hyperbolic 

OPERATORS 

Let L be a normally hyperbolic operator on a globally hyperbolic spacetime X acting on 
the sections of the vector bundle E. Then there exists a unique connection V : C°°{X‘,E) —)■ 
C'°°(A; E 0 T*X) and a potential V € C°°(X; End(F;)) such that 

L = A^ + V, 

where = —Trg(V^®^*^o V) is the connection Laplacian, see [7, Prop. 3.1]. Here 

is the connection on E®T*X induced by the connection on E and the Levi-Civita connection 

on T*X. 

Note that L is of real principal type and the subprincipal symbol sub(L) of L coincides 
in local coordinates with the connection one form of L (see Prop. 3.1 in [17] for a detailed 
explanation). 

Let S C A be a Cauchy hypersurface and let u be its past directed unit-normal vector 
field. Then it is well known that the Cauchy problem for L is well posed and can be solved 
in the sense that there exist continuous operators To : C^{T,;E) — C°°{X;E) and 7i : 
C'“(S; E) —>■ C'°°(A; E) such that for any g, f & C^(S; E) the section u = To/ + Tig is the 
unique solution of the initial value problem 

Lu = 0, ujs =/, {-Tuu)\j]=g. 
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The solution operators 7o,7i are well known to be Fourier integral operators. Unfortu¬ 
nately we were able to find the statement of this fact only for the case of scalar operators. 
Therefore, for the convenience of the reader, we revise the theory in the case of operators 
acting on vector bundles. 

For manifolds V, Z and vector bundles F ^ Z and E ^ Y suppose that C gT*Y x T*Z 
is a homogeneous canonical relation. Let C' = {{y,^,z,rj) \ G C}. Denote by 

jm^Y X Z,C;Hom.{F, E)) the set of Fourier integral operators with canonical relation C. 
These are operators whose Schwartz kernels take values in the bundle E Kl [F* <Z Di(Z)) and 
can be represented locally by an oscillatory integral with a polyhomogeneous matrix valued 
symbol and a scalar phase function that locally parametrizes C. Note that on C we have 
the Keller-Maslov line bundle Me- The principal symbol ajj of a Fourier integral operator 
U G I^(Y X Z, C; Hom(F, E)) is a section in the bundle 

j*c {EMF*)^Mc, 

where jc ■ C ^ Y x Z is the inclusion composed with the projection. A choice of phase 
function that parametrizes the Lagrangian manifold C' gives a local trivialisation of the 
Keller-Maslov line bundle. 

Theorem A.l. Let L as above a normally hyperbolic operator acting on sections of a vector 
bundle E over a globally hyperbolic spacetime. Let Y be a Cauchy hypersurface and let To 
and 7i he the operators solving the Cauchy problem. Then To and Ti are Fourier integral 
operators. More precisely, 

To G I-^X X Y;C]Rom{E,E)) 

Ti G r^X X Y;C;Rom{E,E)), 

where C is the following canonical relation 

C = {{xT,y,r]) G T*X x T*Y \ (x,0 ~ (y,r?)}. 

Here (x,^) ~ {y,y) if and only if there exists a vector fj in T*X such that (x,^) and {y,fj) 
are in the same orbit of the geodesic flow on X and the pull back of {y, fj) to Y coincides with 

{y,v)- 

The statement of the theorem can be implied directly from the calculus developed in [15] 
applied to matrix valued operators as demonstrated by Dencker in [10]. A more direct proof 
for scalar valued operators can be found in [11, Ch. 5.1]. This proof carries over literally 
to the case of operators of real principal type and scalar valued principal symbol. A simple 
computation shows that the characteristic Hamiltonian flow of the principal symbol is the 
geodesic flow in lightlike direction. The subprincipal symbol is the connection one form and 
thus the transport equation (5.1.11) in [11] becomes the differential equation for parallel 
transport along lightlike geodesics. 

Note that in this case there is a canonical choice of trivialisation of the Keller-Maslov line 
bundle obtained by choosing the phase function (j)(x,y,t,t',f) that satisfies (j){x,y,ti,tiT) = 
(x — yT). Assuming the Keller-Maslov bundle has been trivialized this way the principal 
symbols aj-^ and cjt-^ of To and Ti are given as follows: 

a''j^(^xT, y, Tj') F('j^ 

O'Jl {xT, y, rj) (y,r]),(x,^) 1 
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where is the operator of parallel transport along the lightlike geodesic connecting 

(x,^) and {y,ri) by the connection V. 
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